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Abstract 

We study the problem of recovering the initial data of the two dimensional wave equation from values of its 
solution on the boundary dVl of a smooth convex bounded domain VL C R 2 . As a main result we establish 
back-projection type inversion formulas that recover any initial data with support in Q modulo an explicitly 
computed smoothing integral operator JCq. For circular and elliptical domains the operator /Co is shown to 
vanish identically and hence we establish exact inversion formulas of the back-projection type in these cases. 
Similar results are obtained for recovering a function from its mean values over circles with centers on dfl. 
Both reconstruction problems are, amongst others, essential for the hybrid imaging modalities photoacoustic 
and thermoacoustic tomography. 
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1. Introduction 

Let S] C K 2 be a convex bounded domain in the plane with smooth boundary dfl. Suppose that 
/: R 2 — > R is a smooth function that vanishes outside 51, and denote by u: R 2 x (0, oo) — » R the solution of 
the initial value problem 

(<9 t 2 - A x ) u (x, t)=0 for (x, t) e R 2 x (0, oo) 

u(x,0) = f(x) for x € R 2 . (1.1) 

(d t u) (x, 0) = for x e R 2 

In this paper we study the problem of recovering the unknown initial data / from the values of the solution 
u on the observation surface dQ, x (0, oo). Because the solution of (jl.lj) can be expressed in terms of 
circular means (and vice versa), the inversion of the wave equation is basically equivalent to the problem of 
recovering the function / from its averages over circles with centers on dfl. Both reconstruction problems 
are essential for the novel hybrid imaging methods photoacoustic tomography (PAT) and thermoacoustic 
tomography (TAT). The standard setups in PAT/TAT using point-like detectors require the inversion of 
the wave equation in three spatial dimensions [l], Q. The two dimensional version considered in this paper 
arises in a variant of PAT/TAT that uses linear integrating detectors instead of point-like ones (see [3j, |j|). 
Inversion of circular means and the wave equation is also relevant for other imaging modalities, such as 
SONAR Q or ultrasound tomography Q. 

In the recent years, many reconstruction techniques for the wave inversion and the inversion from circular 
means (or spherical means in higher dimensions) have been developed. These techniques can be classified in 
iterative reconstruction methods (see 0, 0, 0, S El ) > model based time reversal (see [h, 12, 1^, H|), Fourier 
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domain algorithms (see [H, [TU [3, 17 , Il8l. [Tit') , and algorithms based on explicit reconstruction formulas of 
the back-projection type (see yl 12, l2dL [2ll I22I 23 , 24, 25, 26, The back-projection approach is 

particularly appealing since it is theoretically exact, stable with respect to data and modeling imperfections, 
mathematically elegant, and quite straightforward to implement numerically. Until recently, however, exact 
back-projection type reconstruction formulas were only known for the cases where dCl is a spherical, planar 
or cylindrical surface (in three spatial dimensions) and a circle or a line (in two spatial dimensions). Many 
researchers even believed that exact reconstruction formulas of the back-projection type exist only in such 
cases. Only very recently the so called universal back-projection formula (originally introduced to PAT/TAT 
in [27} ) has been shown in [25| to provide theoretically exact reconstruction for ellipsoids in R 3 . In the same 
paper it has further be shown that for general convex domains in R 3 the universal back-projection formula 
is exact modulo an explicitly given smoothing integral operator. 

In this paper we establish back-projection type reconstruction formulas for the inversion of the two 
dimensional wave equation on convex domains CI C R 2 . As their counterparts in three spatial dimensions, the 
derived formulas provide exact reconstruction modulo an explicitly computed smoothing integral operator 
/Co. We further show that for elliptical and circular domains the operator JCq vanishes and therefore we 
obtain exact back-projection type reconstruction formulas in these cases. The same type of results are 
derived for reconstructing a function from its mean values over circles with centers on the boundary of a 
convex domain in the plane. Note that our reconstruction formulas as well as the proofs differ from the ones 
given in [25[ . This partially accounts for the fact that the three dimensional wave equation satisfies a strong 
form of Huygens principle, whereas the two dimensional wave equation does not. 

1.1. Notation 

Before presenting our main results we introduce some notation that will be used throughout this article. 

For any k £ N U {00}, we denote by C* (CI) the set of all fc-times continuously differentiable functions 
/ : R 2 — > R which have support in Cl. For given / £ (Cl) we denote by Uf : R 2 x (0, 00) — » R the solution 
of the wave equation (jl.ll) with initial data /. Further, for some integrable function /: R 2 — > R, the circular 
mean transform Aif: R 2 x (0,oo) — > R is defined by 

(Mf)(x,r) := — [ f(x + ruj)duj for (x, r) £ R 2 x (0, 00) . (1.2) 
2tt J s i 

Here and below S 1 :— {2; € R 2 : |x| = l} denotes the unit sphere in R 2 with \x\ being the Euclidian norm 
of 1 e I 2 . The corresponding inner product of two points xo, x\ £ R 2 will be denoted by Xo ■ x\. Similarly, 
the (classical) Radon transform IZf : S 1 x R — > R is defined by 

(Rf)(n,a) := [ f (an + bn^db for (n,a)£S 1 xR, 
Jr 

where n 1 - £ S 1 denotes a unit vector orthogonal to n. The derivative of a function (j>: S 1 x R — >• R with 
respect to the second variable will be denoted by (d a (j)) (n,a), and CH a 4>) ( n >°) will be used to denote the 
Hilbcrt transform of <j> in the second argument, defined as the distributional convolution with l/(7ra). 

We further write xn : R 2 - > R f° r the characteristic function of the domain CI (taking the value one inside 
CI and zero outside) and set 

-I \ ^l-^O -/ \ 1 Nil 2 - No| 2 2 

n(xux ) ■= 1 r, a(xi,x Q ) := - — — for xi ^ x S R . (1.3) 

\xi-x \ 2 laii-afol 

Finally, for any / £ (CI) and any xo £ CI, we define 

lir ,w x 1 f v {dlUaKxn) (n(x 1 ,x ),a(x 1 ,x )) 

(>Cn f) (xo) := — / f(xi)- '—. ; dxi . 1.4 

8tt J n \xi-x \ 

It can be readily verified that the line £ (xi,Xq) — {x £ R 2 : n (x\, xq) ■ x — a (x\,xq)} consists of all points 
having the same distance between xq and x\; compare with Figure fTTTI 



2 



Figure 1.1: Let £(xi,xq) denote the line of all points having the same distance between x\ and xq. Then n(xi,xo) = 
(xi — xq)/\xi — xo\ is the normal vector to £ (x\, xo) and a (x\, xo) = (\xi\ 2 — |^o| 2 )/ (2|xi — xq\) is the oriented distance of 
the line £(xi,xq) from the origin. 



1.2. Main results 

Our first pair of results states that a back-projection type inversion formula applied to the solution of 
the wave equation (jl.ip recovers any initial data / modulo the term /Co/- 

Theorem 1.1 (Inversion of the wave equation). Let CI C R 2 be a convex bounded domain with smooth 
boundary dCl. Then, for every function f £ (Ct) and every reconstruction point xq £ CI, we have 

/W-(MW = -v I0 .f i/« f W){x,t) df , (Uu)- a - } 



TV 



\ x ~ x °\ \Jt 2 ~\X- XqY 



f (x ) - (/Co/) (x Q ) = - / u x ■ (x - x) \ I U ^^ t \ At I As (x) . ( 1 .(i) 



n Jan \ J\x-x \ . Uz _ u _ t„i 2 



\x-x a \ ^j t 2 -\ x - x \ 

Here Uf is the solution of /Co/ is defined by ^l-4\ l, df is the usual arc length measure on dCl, v x 

denotes the outwards pointing unit normal to dCl, and V^- denotes the divergence with respect to xq. 

Proof. See Section El □ 

We also prove the following corresponding formulas for recovering a function from its mean values over 
circles centered on the boundary dCl. 

Theorem 1.2 (Inversion from circular means). Let CI C K 2 be a convex bounded domain with smooth 
boundary dCl. Then, for every f £ (CI) and every xo £ CI, we have 

f (xo) - (/Co/) (xo) = r V. D • / u x ( r r ! Mf) {X,r l d M ds (*) . (1-7) 



Jen \Jo r 2 — \x — xq\ 



/(•'•„)- lk\,_. /■)(./•„! = ~ / ^ ■ (xo - x) [ r {d J Mf) ^ T } 2 dr ) ds (g) . (1.8) 



^ -/ao \Jo r 2 — |x — x \ 



Here both inner integrals have to be taken in the principal value sense, Aif is the circular mean transform 
defined by \1.2\) . and /Co 7 v x , d^ and V Xo are as in Theorem ] 

Proof. See Section |4] □ 



Remark 1.3 (Smoothing effect of ICn). // O C R 2 is smooth and convex, then the Radon transform of 
Xn is smooth except for such pairs of (n, a) £ S 1 x R where the line {x £ R 2 : n ■ x = a} is tangential to 
the boundary dfl. Notice further, that for x\ 7^ xq € £1 the line t{x\,Xo) is never tangential to dtt; see 
Figure ["PI Since the operator d^H a preserves the locations of the singularities oflZxn this shows that the 
kernel in |1.^| ) has at most a weak singularity proportional to 1/ \x± — Xq\ and hence is at least smoothing 
by one degree. 

For speciai domains, the integral operator ICq may vanish, in which cases Theorems 11.11 and 11.21 provide 
exact reconstruction formulas for the inversion of the wave equation and the inversion from circular means, 
respectively. We verify that this is indeed the case for circular and elliptical domains. 

Theorem 1.4 (Exact inversion for circular and elliptical domains). Let £1 C R 2 be a circular or elliptical 
domain and let f: R 2 — > R be a C°° function with support in fl. Then /Cn/ vanishes identically on Q. In 
particular, the following hold: 
(a) The function f can be recovered from the solution of the wave equation M.l\) by means of either of the 
following formulas: 



f(xo) = - V Xo ■ [ Vx \ [ 

K JdQ \J\x-xo 



{Uf)M :d*ld*(x), (1.9) 



I yjt*-\x-x f 



f(x )=- [ Vx -{x ~x)[r ^MMdt| ( u(,). (i.io) 

11 Jm \J\x-x„\ ^J t 2 _ \ X - Xq \ 2 

(b) The function f can be recovered from the circular averages Mf defined in by means of either of 

the following formulas: 



f (x ) = - V X0 ■ / v x \ l 
n Jan \Jo 



r(Mf) (x,r) , \ , , . 

— ^ — JJV ' 2 dr ds (x) , (1.11) 
r 2 -\x-x \ I 



f (a*) = - [ v*-(x -x)(r {d : M p {x ' T l dr) ds (x) . (1.12) 
T Jan \Jo r 2 - \x - x \ J 

Here both inner integrals have to be taken in the principal value sense. 
Proof. See Section [S] □ 

1.3. Prior work and innovations 

In the case that is a spherical domain in d > 2 dimension, various exact back-projection type inversion 
formulas for the inversion of the wave equation and the inversion from spherical means have been derived 
in [3, 12, 21, 23, 27 j. In particular, the inversion formula Q coincides with our formula (jl.101) for the wave 



inversion, and the formula of [23j for the inversion from circular means (if rewritten as in [28| ; see [29j for a 
different derivation) coincides with our formula (|1.11[) . However, in [3I. [23j[ these formulas are neither shown 
be exact for elliptical domains, nor are they investigated for more general domains. 

Inversion formulas for the spherical mean transform on elliptical domains have been derived recently 
in [2^|. The methods and results there are different from ours and no statements are made for general 
convex domains. Even the reconstruction formula of j2(| for ellipses differs from each of our formulas (jl.ll[> . 
(|1.12[) . Our results are more closely related to the one of [25[, where corresponding results have derived for 
convex domains in R 3 . Note further that in three spatial dimensions a statement similar to our Theorem ll.il 
is also present in (27 1 (w here, however, no explicit expression for the operator /Co has been derived). Anyway, 
none of the papers 



25l . |27| considers the two-dimensional case. 
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1.4-. Outline 

The following sections are devoted to the proofs of the theorems presented in Section 11.21 In Section [5] 
we derive auxiliary statements that we require for the proofs of these results. In Section [3] we then derive 
the formulas for the wave inversion claimed in Theorem 11.11 and in Section [4] we establish the corresponding 
results for the inversion from circular means. The cases of circular and elliptical domains are considered in 
Section where we show that the operator ICq vanishes in these cases and hence we establish the exact 
inversion formulas claimed in Theorem II .41 The paper concludes with a short discussion in Section [5] 



2. Auxiliary results 



For the following considerations recall the definitions of the circular mean transform A4, the Radon 
transform 1Z, and the operator U that maps any initial data / € (fi) to the solution of (|1.1|) ; see 
Section 11.11 We will further make use of the operator V, which maps any / 6 (Q) to the solution 



x (0, oo) — ► R of the initial value problem 



' (<9 2 - A x ) v (x, t) = for (x, t) € R 2 x (0, oo) 



u(a:,0) = 
(d t v) (x,0) = f(x) 



for i£l 

for x em 2 



(2.1) 



One easily recognizes that U and V are related by the identity Uf — dtVf. Indeed, if v solves (|2.1[) . then 
by definition u = dtv satisfies the wave equation (9 2 — A x )u = and the initial condition u (x, 0) = / (0). 
Further, the wave equation yields (d t u) (x, 0) = (d^v) (x, 0) = (A x v) (x, 0) = which shows that u is a 
solution of (jl.ip . 

Next we proof a simple integral identity that we require for the following considerations. 



Lemma 2.1. Suppose that f <E (f2) and that . 
outside Cl. Then, for every x\ € fl, we have 



4>{x) 



(d ro Mf) (x,r )ln 



1 

2^ 



l I u 

Tq-\X-Xx\ 

/Oo) 



\X\ - X \ 

Here ft = h (xi, xq) and a = a (x\, xq) are defined by Equation fl.3\) 



R is continuously differentiable on Q and vanishes 
<!/'o J da; 

{d a 1Z4>) (ft, a) In |2 \x\ — Xq\ (a — a) \ da I dxg • (2.2) 



Proof. We first verify (|2.2|) with $ ( ■ ) in place of In | • | , where $ : R — >• 
integrable. Hence we show that for every x\ £ il the following holds: 



is continuously differentiable and 



(x) (J (d ro Mf) 0, n>) * (rg — |ar— £i| 2 ) dr ) 

/Oo) 



d.r 



1 
2^ 



|xi - ar | 



(9 o ft0) (ft, a) $ (2 |xi - a; | (a - a)) da dx . (2.3) 



After performing one integration by parts, using the definition of the circular mean transform, and intro- 
ducing polar coordinates around the center x afterwards, the inner integral on the left hand side of (|2.3I) 
can be written as 



(d ro Mf) (x, r„) $ (rg -\x- Z!| 2 ) dr = -2 jf r (A*/) (a, r ) (rg - |x - xx| 2 ) dr 

rof (x + r u ) $' (r 2 , - | a; - xi| 2 ) dw dr 



- / / Oo) (|& - x | 2 - |.t - £i| 2 ) dx 
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By straight forward computation one verifies that \x — xo| 2 — |x — xi\ 2 = 2 (xi — xq) ■ [x — (xi + xq) /2). 
Therefore, one application of Fubini's theorem, the definition of the Radon transform and the definitions of 
h = n (xi , xq) and a = a (xi , xq) yield 



f 4>{x) r (d ra Mf)^,r )<f(r 2 ~\x^x 1 \ 2 )dr 

R 2 JO K 7 



d.r 



— If ipo) ( <t> ( x ) ( 2 {x\ - x ) ■ {x - (xi - x ) /2)) dxdx 
w Jn Jr 2 

-— J f(xo) (^J J ^(ah + bh 1 -) $' (2 \x\ - ar | (a - a)) dbdaj dx 
f (x Q ) [ / (Tty) (ft, a) $' (2 |xj - x | (a - a)) da j dx 



Finally, after performing one further integration by parts in the inner integral, the last expression is seen to 
be equal to 

1 f /Oo) ( f (fl a 7ty)(a j0 )$( 2 |a;i - x | (a - a)) da ] dx . 



2tt Fi ~ 

This shows that (|2.3p indeed holds for all f eC 1 (M) n L 1 (R). In order to verify the corresponding identity 
for In | • | in place of $ we may write In | • | as the L 1 -limit of continuously differentiable functions $„. 
Application of (|2.3p with $„ in place of $ and taking the limit n — > oo afterwards then implies the desired 
identity (|2~2]l , □ 



The next Lemma is the key for the results in this paper. 
Lemma 2.2. For all f,g 6 (fl) we have 



(Uf)(x,t) (Vg) (x,t) dt)dx 



= -^~ / /WffOi) i j dx dxi. 2.4 

&nJ n J n |xi-x | 

.ffere "H a «s the Hilbert transform in the second argument, \a * s the characteristic function of Q, and n = 
n(xi,xo) and a = a(xi,Xo) are defined by Equation 11. 3\) . 

Proof. First recall that the explicit expressions for the solutions of the initial value problems (|2.ip and 
(see, for example, [3(J p. 134]) are given by 

s / ,s f l r (Mg) (x,r) , 

„, fU +A fl /"* r(Mf) (x,r) ^ 
(Uf)(x,t) = d t J Q -j== dr. 

After performing one integration by parts and differentiating under the integral sign, the expression for Uf 
can be rewritten as 



(///Mr./! = u, j (d r Mf)(x,r)Vt 2 -r-dr= j (d r Mf) (x,r) J—^ dr. (2.5) 



G 



These identities and Fubini's theorem yield 



(Uf) (x,t) (Vg) (x,t)dt 



T-i-oo 



= lim / (riMg) (x, ri 



T ^°°J0 JO 

The inner integral in the above expression evaluates to 



<{r ,n} \Jt 2 - rl^Jt 2 - r\ 



In ( JT* -rl + JT* -rA-lte (|r 2 - r 2 |) . (2.6) 



Hence we obtain 

Id* 



/•OO 

/ (Uf)(x,t)(Vg)(x,t), 
Jo 

(d ro Mf) (x, r ) (riXg) (x, n) In (jT*-rl + ,Jt*-j^ dr dn 



poo poo 

lim 

T-^oo 



OO /-OO 

/ (d ro Mf) (x,r ) (nMg) (x,ri)\n\r% - rf\ dr drj 
o Jo 



After one integration by parts and performing the limit T — > oo afterwards, the first term on the 
right hand side is seen to vanish. Therefore, by integrating over the variable x, interchanging the order of 
integration and introducing polar coordinates, we obtain 



Q JO 



g{xi)\ / Xn{x) / (d ro Mf) (x,r )ln r 2 - \x - x x [ 



(Uf) (x,t) (Vg) (x,t) dtdx 
1 

477 

According to Lemma |2~T1 the double integral in brackets equals 

/(so) 



drodx I dxi 



2tt J n \xi - x \ 



(daJtxn) (ft, o) In 1 2 |xi - x \ (a - a)\ da dx ■ 



Next note that the distributional derivative of ln|o| is P.V.^. After recalling the definition of the Hilbert 
transform we therefore obtain 



n Jo 



(Uf) (x,t) (Vg)(x,t) dtdx 

f(xo)g(xi) 



(daTi-Xfi) (ft, a) In \2 |xi - x \ (a — a)\da) dx dx 



1 f [ f( x o)g(xi) f f, w .da , , 

P.V. / (Uxn) (n,a) )dx dxi 



8tt 2 J q J q \xi-Xq 
1 f f I(xa)g(xi 



(HaTi-Xn) (ft, a) dx dxi 



Hence we have verified (|2.4I) which concludes the proof □ 
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We are now ready to derive the following auxiliary theorem from which we will extract all formulas 
presented in the introduction. 

Theorem 2.3. For all f,g€C™ (fi) we have 

»x ■ ( I (Uf) (x, t) (V x Vg) (x, t) At) As (x) 



on 



= / / (xo) 9 (xq) dx - / (K,nf)(xo)g(x )Ax . (2.7) 



Here U and V denote the solution operators of the initial value problems and \2. 1)) , respectively, K,nf 

is defined by and v x denotes the outward pointing unit normal to d£l. 

Proof. Application of Greens second identity to the functions (Uf) ( • , t) and (Vg) ( ■ , t) with fixed t followed 
by an integration over the temporal variable yields 



((Vg) (x, t) (A x Uf) (x, t) - (Uf) (x, t) (A x Vg) (x, t)) Ax J At 

v x ((Vg) (x, t) (V x Uf) (x, t) - (Uf) (x, t) (V x Vg) (x, t)) As (x) ) At . 



ion 

Since V f and Ug are both solutions of the wave equation we have the equalities A x Uf = dfUf and A x Vg = 
dfVg. Together with the relation (Vg) (V x Uf) — V ' x (VgUf) — (V x Vg) (Uf) and one application of the 
divergence theorem this further implies 



n wo 



((Vg) (x, t) (d 2 t Uf) (x, t) - (Uf) (x, t) (d 2 t Vg) (x, t)) At Ax 



A x ((Vg) (x, t) (Uf) (x, t)) Ax) At -2 [ v x ((Uf) (x, t) (V x Vg) (x, t)) As (x) At . 

J Jo \Jon J 

After integrating the inner integrals on the right hand side by parts and using the initial conditions 
(Uf)(-,0) = f, (dtUf) (■ ,0) = and (d t Vg)(-,0) = g, the expression on the left hand side evaluates 
to J n f(x )g (x ) Ax which yields 

f(xo)g(x )Ax = f I A x ((Uf)(x,t)(Vg)(x,t))AtAx 
Jn Jo 



v x \ (Uf) (x, t) (V x Vg) (x, t) AtAs (x) . (2.8) 
on Jo 



Next notice that 



A, ((Uf) (Vg)) = (A X U f) (Vg) + 2 (V x Vf) ■ (V x Ug) + (Uf) (A x Vg) 
= (UA x f) (Vg) + 2 (VV x f) ■ (UV x g) + (Uf) (VA x g) , 

where VV X / (and similarly UV ' x g) is defined by component-wise application of V to V ' x f '. Hence in view of 
Lemma 12.21 and by repeated application of the divergence theorem we obtain 

A x ((Uf) (x,t) (Vg) (x,t))AtAx 

Iff, .2,,, x , ^ (n a nxu)(n(x 1 ,x ) ,a(x 1 ,x )) 

= ~^~ / Vxo+VxJ / (X )g Xi ; ; Ax Ax! 

^ JnJn |xi-x | 

1 f f t ( \ i uu u \ 2 f(rla'R-xn)(n(x 1 ,xo) ,a(x u x ))\ 

= ~^~ / f( x o)9(xi)(Vx +V S1 ) ; ; dx dxi. 

8tt JnJn \ xi - x J 



n Jo 



As easily verified, the following relations hold true: 

(V Xl +'V Xo )n(x 1 ,x ) = 0, 

(V Xl +V x „)a(x 1 ,x ) = (xi - x ) / \xi -x \ , 

(V 2;i +V Xo )|xi-a;or 1 =0, 
(V Xl +V ;Eo )-(^i-a;o)=0. 



n Jo 



These relations imply 

(V Xl + V X0 ) — — — = (V X1 + V Xo ) — -2 (d a n a nxn) [n, a) = — : , 

8n\x 1 -x Q \ 87r|xi-a;Q| 8?r - x | 

and therefore 

A x ((Uf) (x,t) (Vg) (x,t))dtdx 

1 f f tt \ I \ ( d l^g' R 'Xn)(n(x 1 ,xo),a(x 1 ,x j) f 

= / / j{ x o)g(xi) j j dx dxi= / {K.nf)(x 1 )g(xi)dx 1 . 

^JnJn \xi-x \ Jq 

Inserting the last displayed equation into Equation (|2.8|l hnally yields the claimed identity (|2.7[) . □ 

3. Proof of Theorem 11.11 

In this section we establish the reconstruction formulas stated in Theorem 11.11 as a consequence of 
Theorem 12.31 and the following Lemma that is easy to establish. 



Lemma 3.1. For every f,g£ (Q) we have 



v x - / {Uf){x,t)(V x Vg)(x,t)dtds(x) 
an Jo 



V Xg ■ / v x / = dtds (x) g (xq) dx 



2- 



an J\x-x \ Jt 2 - |ai — ac Q | 



Proof. The identity V x Vg = W x g (again W x g is defined component-wise) and the explicit expression for 
the solution of the wave equation (|2.ip imply 

/>oo 

v x - \ (Uf)(x,t){V x Vg)(x,t)dtds(x) 
an Jo 

1 / v x ■ r (Uf) (x, t) [ (V Xo g) (x a ) X j lx ~ Xol<t} dx a dtds (x) 



2n 



an Jo 



\Jt 2 - \x - x Y 



= 7T I (Vx„s)(a:o)- / "* T mM dtds(x)dx 

= -T-f ( V *°7 "* I" dtds(x))g(x )dx Q . 
27T JW\ Jon J\x-x \ ^ t 2_\ x _ XQ f J 

This is already the desired identity and concludes the proof of Lemma 13.11 □ 
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3.1. Proof of formula \1.5}) 

By Theorem 12.31 and Lemma T3.1[ for all /, g G (f2) the following identity holds true: 



/ (x ) g (x ) dx - / (/Co/)(x )sf(xo)dx . 



7T ./iB2 



This implies that (|1.5[) holds almost everywhere. Since both sides of (|1.5p are continuous functions the 
equality must also hold pointwise. 

3.2. Proof of formula [FQ) 

The second identity given in Theorem 11.11 is an easy corollary of the first one just established. Indeed, 
from the chain rule and one integration by parts we obtain 



7T 



on 



x-x \ ^J t 2 -\ x - Xq\' 
1. 



-V^ • / v. 

7T 



<9S1 */ 1 ic — 



l. 



= V X0 • / V. 

7T 



BVl J\x— xq 



d t (^ 2 -|x-x | 2 ) (J -1 !*/) (x,t) dtds(x) 
^2- |x-x | 2 (dtt^Uf) (x, t) dt ds (x) 



/ ^ x • (xq — x) ( f - ' = = ' = dt I ds (x) 



7T 



^-^l ^ 2 - |x - X | : 



In view of fj 1 . 5f) . this verifies formula (|1.6p . 

4. Proof of Theorem [TT2] 

In this section we derive the inversion formulas for the circular mean transform stated in Theorem 11.21 
The proofs will be based on the inversion formula (| 1 . 5|) for the inversion of the wave equation and the explicit 
expression (I2.5[) for the solution of (jl.ip in terms of the circular means A4f. 

4-1. Proof of formula Ji.Tj ) 

The inversion formula (|1.5[) for the wave equation and the explicit expression (|2.5p for the solution of (|l.ip 
imply that 



-V, -/ «J hm / T f ] d, ( x) . (4.1) 
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After changing the order of integration and using Equation (|2.6[) . the inner double integral evaluates to 



oo pT 



lun I i), (Mf)(.r. r) / ^L= dr 







(9 rJ M/)(ac,r)ln r 2 - |x-x r dr 



Together with Equation (|4.1I) and one integration by parts (using that the distributional derivative of In \ r\ 
is 1/r) this yields 

f(xo)-(JCnf)(x ) = -^-V Xo -[ v x ( (d r Mf)(x,r)\n\r 2 - \x - x \ 2 \ drds(x) 
27T Jon Jo 

= iv,,/ JriMMXi,), ,,2, 

^ ./an \yo — |x — xo| y 

where the last integral is taken in the principal value sense. The last identity obviously coincides with the 
inversion formula (|1.7[) . 



4.2. Proof of formula U.8}) 

The second inversion formula (|1.8[) for the circular mean transform M could be obtained from (|1.6[) in 
way similar (but slightly more involved) to the derivation of (jl.7p presented above. However, it is simpler 
to derive formula (|1.8p directly from (|1.7p . 

One integration by parts and interchanging the order of differentiating and integration yields 

l -V X0 .(( Vx r rMf ^\ 2 dr)d S (x) 
n \Jdn Jo r 2 ~\x-x \ J 

= -^V Xn - f v x \ {d r Mf){x,r)\n\r 2 -\x-x \ 2 \drds{x) 
27r Jon Jo 



= -^T / v *- I V So (lii|r 2 -|a:-a:o| 2 |)(fli.-M/)(a:,r)drd«(a:) 
ZlT Jan Jo 



n JdQ, \Jo r 2 — \x — x 0l / 

In view of Q1.7P this verifies the formula (|1.8p . 

5. Proof of Theorem 11.41 

In this section we verify the exact reconstruction formulas presented in Thcorcm ll.4l in the case that the 
domain 57 is a disc or an elliptical domain. According to Theorems 11.11 and 11.21 it is sufficient to show that 

{dlUaRxn) (n(x-L,xo) ,a(x-L,x )) =0 for allxi,x G Q, 

where n (x\, xq) and d {x%, xq) are defined by (|1.3p . 

5.1. Circular domains 

We first consider the special case where £1 = D is a disc in the plane and we assume without loss of 
generality that D = {x £ R 2 : \x\ = l} is the unit disc centered at the origin. 
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The Radon transform of \d equals 2y/l — a 2 x{|a| 2 < l}. The Hilbert transform of the function 
y/l — a 2 x{ \ a \ 2 < l} is known (see, for example, [3l|, Table 13.11]) and yields 



-a — y/ a 2 — 1 if a < —1 
Vl{„\,,)l».a} : 2 { -a if - Ks < 1 . 

-a + V.5 2 - 1 if a > 1 

This in particular implies that {d 2/ H a XB) (n,a) = for \a\ < 1. Because we have a(xi,xo) < 1 for all 
xi,xq G D, this implies Theorem II .41 for the case of circular domains. 

5.2. Elliptical domains 

Now suppose that f2 = E is an elliptical domain. We may assume without loss of generality that 

E = j(x, y) e M 2 : x 2 + |^ < 1 j for some b > . 

We then have Xn(x, y) = Xd {x,y/b), where D is the unit disc as above. For any integrable func- 
tion ip: E 2 — >• K and any invertible matrix B £ R 2x2 one can easily verify the identity (7Zip B ) in, a) = 

Tgr^T (^V 7 ) (n/|i? T n.|, a/|S T n|) with (x) :— ip (B~ 1 x). Consequently, after writing 9(a) — (cos a, sin a) 
we obtain 

6(a) g 



bn XD 

(n XE )(0(a),a) = 



\J cos 2 a+b 2 sin 2 a ' \J cos 2 a+b 2 sin 2 a 



\/ cos 2 a + 6 2 sin 2 
and therefore 

bidtUaTlXD) 1 ~ 

(d 2 a n a Kx E ) (0 (a) , a) = 



y cos 2 sin 2 a ' ^ cos 2 a+b 2 sin 2 a 



cos 2 a + 6 2 sin 2 a) 3 ^ 2 



For a?o,Xi € £1 write n(xi,xo) = 0(a). Then we have a(xi,xo) /y cos 2 a + a 2 sin 2 d < 1. From Subsec- 
tion l5.ll we already know that (9 2 H q 7?.x_d) (n,a) = for every \a\ < 1, and therefore obtain 

(djUaTi-XE) (n (xi, x ) , a (xi,x )) = for all x 7^ xi e O , 
which establishes Theorem II .41 for the general case of elliptic domains. 



6. Conclusion 

In this paper we derived inversion formulas of the back-projection type that recover the initial data of 
the wave equation from its solution on the boundary dfl of a convex domain modulo the smooth term 

fir , w \ 1 f . {dlUaTLxn) (n(x 1 ,x ) ,a(xi,x )) 

(K(lf) (X ) = — / / (Xl) ^ '- : dXl 

8tt J n |xi-x | 
(see Theorem II. ip . In the case of circular and elliptical domains the operator /Co has been shown to vanish 



identically which yields to exact reconstruction formulas (see Theorem ll.4l Item (a) I in these cases. Corre- 
sponding statements have been derived for the inversion of the circular mean transform (see Theorem 11.21 
and Theorem ITU Itcm|(b)[). 

We note that all formulas derived in this paper can be implemented in a quite straight forward manner 
following the derivations in 0, [2l[ . We do not give details here and refer the interested reader to the papers 
[3H2III for detailed derivations of discrete back-projection type algorithms. A numerical reconstruction based 
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Figure 6.1: Numerical reconstruction using inversion formula 11. 911 on the elliptical domain E = {(x,y) : x 2 + (y/0.8) 2 < 1}. 
Left: Initial data /. Middle: Simulated data Uf. Right: Numerical reconstruction using inversion formula 11. 9H . 



on the inversion formula (|1.9p on the elliptical domain E = {(x, y) : x 2 + (y/0.8) 2 < 1} is shown in Figure RTT1 
It can be seen that except for some smoothing effects at boundaries (due to the numerical implementation), 
the initial data is recovered almost perfectly. 

With G denoting the diverging fundamental solution of the two dimensional wave equation and with u 
denoting the solution of (jl.ip , our inversion operator in (11.51) and (11. 9| can be written in the form 

/ubp (x ) = 2 V I0 • / v x l G(x - x ,t)u(x,t)dtds(x) . 
Jan Jo 

The analog of this expression in three spatial dimensions (where G and u are replaced by the three dimen- 
sional fundamental solution and the solution of the three dimensional wave equation, respectively) is the 
so-called universal back-projection introduced in the context of photoacoustic tomography in [27j . In this 
paper it has also been shown that the universal back-projection formula exactly recovers the initial data of 
the three dimensional wave equation on spherical domains. This result has been extended to ellipsoids in 
M 3 in 25]. In the latter paper, it is further shown that for any smooth convex bounded domain Q C M 3 , the 
universal back-projection formula recovers the initial data of the three dimensional wave equation modulo 
the term 



/V(3D)A/ x 1 f tl s (daKXn) (n(x 1 ,x ),a{x 1 ,x )) 

V ) 16tH J q \xx-xol 2 



dxi 



The identities derived in the current paper demonstrate that results similar to the ones of |25| also holds in 
two spatial dimensions. After initial submission of the present manuscript, these results have actually been 



generalized to arbitrary dimension (see 32]). 
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